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1. Introduction 

Quantum graphs of one-dimensional wires connected at nodes were introduced already 
more than half a century ago to model physical systems. Depending on the envisaged 
application the precise formulation of the models can be quite diverse and ranges from 
solid-state applications to mathematical physics [Tll2llSllllElinilZllSlinilIIIlini[I2!. 
Lately, quantum graphs attracted also the interest of the quantum chaos community 
because they can be viewed as typical and yet relatively simple examples for the large 
class of systems in which classically chaotic dynamics implies universal correlations in the 
semiclassical limit [TnillllinilinillllllHllISlEnilZIlEZI. Up to now we have only a limited 
understanding of the reasons for this universality, and quantum graph models provide a 
valuable opportunity for mathematically rigorous investigations of the phenomenon. In 
particular, for quantum graphs an exact trace formula exists [121 HBl Hlj which is based 
on the periodic orbits of a mixing classical dynamical system. Moreover, it is possible 
to express two-point spectral correlation functions in terms of purely combinatorial 
problems [13 CZl UHl UH 120] and to use this in order to investigate the origin of the 
connection between Random Matrix Theory (RMT) and the underlying classical chaotic 
dynamics. 

By attaching infinite leads at the vertices, we get non-compact graphs, for which a 
scattering theory can be developed |211 12S] • The quantum scattering matrix for such 
systems can be written explicitly, and the corresponding resonances can be calculated 
easily. It is the purpose of this paper to review their statistical properties and analyze 
their dependence on the boundary conditions that we impose at the vertices [201 123 12H1 ■ 
At the same time we apply the reaction-matrix formalism and for the first time we give 
a semiclassical interpretation of this theory in terms of the classical trajectories of the 
graph. 

The paper is structured in the following way. In section |2l the mathematical model 
is introduced and the main definitions are given. Section [SI is devoted to the derivation 
of the scattering matrix for graphs. In section H] we discuss the relation between the S- 
matrix of an open graph and the spectrum of the corresponding closed system (reaction- 
matrix formalism). The resonance width distribution for various types of graphs are 
presented and discussed in section [HI and deviations from RMT predictions [201 1211 12H] 
are analyzed. In the same section we also exhibit and analyze the phenomenon of 
resonance trapping occurring in specific graphs. Our conclusions are summarized in 
section |01 

2. Quantum Graphs: Definitions 

We start by considering a compact graph Q. It consists of V vertices connected by B 
bonds. The number of bonds which emanate from each vertex i defines the valency Vi 
of the corresponding vertex (for simplicity we will allow only a single bond between 
any two vertices). The graph is called v — regular if all the vertices have the same 
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valency v. The total number of bonds is S = | Yh=i "^i- Associated to every graph is its 
connectivity matrix C . It is a square matrix of size V whose matrix elements Cij take 
the values 1 if the vertices i,j are connected with a bond, or otherwise. The bond 
connecting the vertices i and j is denoted by 6 = (i, j), and we use the convention that 
i < j. It will be sometimes convenient to use the "time reversed" notation, where the 
first index is the larger, and b = (j, i) with j > i. We shall also use the directed bonds 
representation, in which b and b are distinguished as two directed bonds conjugated by 
time-reversal. We associate the natural metric to the bonds, so that Xij (xj^i) measures 
the distance from the vertex i (j) along the bond. The length of the bonds are denoted 
by Lb and unless stated otherwise we shall assume that they are rationally independent. 
The mean length is defined by (L) = {l/B)J2b=i^b and in all numerical calculations 
below it will be taken to be 1. In the directed- 6on(i notation Lj, — Ly. Finally, we define 
the total length of the graph C = Y.b=i Lb- 

The scattering graph Q is obtained by adding leads which extend from M(< V) 
vertices to infinity. For simplicity we connect at most one lead to any vertex. The 
valency of these vertices increases to Vi — Vi + 1. The M leads are denoted by the index 
i of the vertex to which they are attached while Xi now measures the distance from the 
vertex along the lead i. 

The Schrodinger operator (with h — 2m = 1) is defined on the graph Q in the 
following way: On the bonds b, the components of the total wave function ^ are 
solutions of the one - dimensional equation 

(-i^-^) <ilb{x) = k'^bix), b^{t,j) (1) 

where Ab (with ^e{Ab) ^ and Ab = —A-^) is a "magnetic vector potential" which breaks 
the time reversal symmetry. In most applications we shall assume that all the AbS are 
equal and the bond index will be dropped. The components of the wave functions on 
the leads, 'ifi{x), are solutions of 

-^^i{x)^e^i{x), i = l,...,M. (2) 

At the vertices, the wavefunction satisfies boundary conditions which ensure current 
conservation. To implement the boundary conditions, the components of the wave 
function on each of the bonds b and the leads i are expressed in terms of counter 
propagating waves with a wave- vector k: 

On the bonds : ^b = abe'^''+^'>^''' + Cbe'^-'^+^b)xb 

On the leads : *i = /^e-^^^^* + O^e^^^^* . (3) 
The amplitudes a^, Cb on the bonds and /j, on the lead are related by 
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These equalities impose the boundary conditions at the vertices. The vertex scattering 
matrices sj*^-,, are Vi x Vi unitary symmetric matrices, and go over all the Vi bonds 
and the lead which emanate from i. The unitarity of S^*) guarantees current conservation 
at each vertex. 

On the right hand side of (0}, the vertex scattering matrix S'^*) was written 
explicitly in terms of the vertex reflection amplitude p*-*-*, the lead-bond transmission 
amplitudes {'rj*''}, and the fj x Vi bond-bond transition matrix which is sub unitary 
(I det o"*^*^ I < 1), due to the coupling to the leads. Vertices which are not coupled to leads 
have p'-*-' = 1, rj*'' = 0, while the bond-bond transition matrix crj*j, is unitary. 

Graphs for which there are no further requirements on the S*^*) shall be referred to 
as generic. In this case the vertex scattering matrices E*^*) are chosen from the ensemble 
of unitary symmetric Vi x Vi random matrices, as explained in ^21- However, often it is 
more convenient to compute the vertex scattering matrices from the requirement that 
the wave function is continuous and satisfies current conservation at all vertices P^ ll4j. 
In this case the resulting E*^*) matrices read fl^ : 

^,,^(i+e-)_ ,(.._(i±£r^^ /»^ (^ + 1(5) 

J,J y J'J J y y 

where Ui = 2arctan;i^. The parameters < Aj < oo are free parameters which 
determine the boundary conditions. We shall refer to the A, as the vertex scattering 
potential. The special case of zero Aj corresponds to Neumann boundary conditions. 
Dirichlet boundary conditions result from A^ = oo. A finite value of Aj introduces a new 
length scale. It is natural therefore, to interpret it in physical terms as a representation 
of a local impurity or an external fields We finally note that the above model can 
be considered as a generalization of the Kronig-Penney model to a multiply connected, 
yet one dimensional manifold. 



3. The S-matrix for Quantum Graphs 

It is convenient to discuss first graphs with leads connected to all the vertices M = 
V. The generalization to an arbitrarily number M < V of leads (channels) is 
straightforward and will be presented at the end of this section. 

To derive the scattering matrix, we first write the bond wave functions using the 
two representations which are conjugated by "time reversal": 

^^bixb) = abe'^''+^''^'"' + Cbe'^-^^^^^'"' = 

Hence, 

= age^^'^+^^^^S afc = qe*(-'=+^^.)^''. (7) 

In other words, but for a phase factor, the outgoing wave from the vertex i in the 
direction j is identical to the incoming wave at j coming from i. 
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Substituting aj, from Eq. ((7j) in Eq. (jD), and solving for Cjj- we get 

r,s ^ ' ''^ ' 

=p«/, + ^r«c,,, (8) 
i' 

where 1 is the 2B x 2B unit matrix. Here, the "bond scattering matrix" Sh is a sub- 
unitary matrix in the 2B dimensional space of directed bonds which propagates the 
wavef unctions. It is defined as Ssik^A) = D{k;A)R, with 

A,y,'(A;, A) = 5,,,,5,-,,e*'=^-+^^'.^^- (9) 

Rji,nm ^n,iCj^iCi^rn(^ji^im' 

D{k,A) is a diagonal unitary matrix which depends only on the metric properties of 
the graph, and provides a phase which is due to free propagation on the bonds. The 
sub-unitary matrix R depends on the connectivity and on the bond-bond transition 
matrices a. It assigns a scattering amplitude for transitions between connected directed 
bonds. R is sub-unitary, since 

V 

Idet^l = n |deta(^)| < 1. (10) 

i=l 

Replacing Ciji in the second of Eqs. (jH)) we get the following relation between the 
outgoing and incoming amplitudes Oi and Ij on the leads: 

= p«/. + E (l - Snik; A))~l^^^^^^^ D^,,,rl^h, . (11) 
Combining (fTTjl for all leads z = 1, . . . , l^, we obtain the unitary V xV scattering matrix 

= k^P^'^ + E ^i^^ (l - Ssik; A)) -'^^^^^^^ D(,,)ri^-). (12) 
Finally, a multiple scattering expansion for the S*— matrix is obtained by substituting 

oo 

{l-SB{k;A))-' = Y.SUk;A) (13) 

n=0 

into Eq. ()12|1 . The expansion (fT!?|l converges absolutely on the real fc— axis since for 
generic graphs the spectrum of S is strictly inside the unit circle (or else, for k + ie, with 
e > arbitrarily small). 

The scattering matrix can be decomposed into two parts S^'^\k) = {S^'^^)k + S-I'''{k) 
which are associated with two well separated time scales of the scattering process. 

can be evaluated by averaging term by term the 5*— matrix resulting from the 
substitution of the multiple scattering expansion (fT^ in (fT^ . Since all terms, but for 
the prompt reflections, are oscillatory functions of k, they average to zero (we assume 
that the fc— window over which the averaging is performed is sufficiently large i.e. larger 
than 27r/Lmin) and thus we get that {S^^^)k = = 5ijp^^\ The fluctuating component 
of the S matrix, S-^''{k), starts by a transmission from the incoming lead i to the bonds 
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(i,r) with transmission amphtudes Tj;^\ The wave gains a phase e*^'^"'""^'')^'' for each 
bond it traverses and a scattering amphtude 5"^*] at each vertex. Eventually the wave is 
transmitted from the bond {s,j) to the lead j with an amplitude t^^\ Explicitly, 

5(y = 5,,,p«+ i3ie^('^''+®*) (14) 

where T^^j is the set of the trajectories on Q which lead from i to j. Bt is the amplitude 
corresponding to a path t and is given by the product of all vertex-scattering matrix 
elements encountered along t. Length and directed length of the path are k = J2bet Lb 
and Qt = J2b& Lb^b, respectively. Thus the scattering amplitude S'f J-* is a sum of a 
large number of partial amplitudes, whose complex interference brings about the typical 
irregular fluctuations of |5'j-^''p function of k. 

We finally comment that the formalism above can be easily modified for graphs 
where not all the vertices are attached to leads. If the vertex / is not attached, one has 
to set = in the definition of S'^'^ The dimension of the scattering matrix 

is then changed accordingly. 



4. Relations between spectrum and S-matrix 
4.1. Secular equation 

There exists an intimate link between the scattering matrix and the spectrum of the 
corresponding closed graph. It manifests the exterior- interior duality for graphs [T^I23j . 
The spectrum of the closed graph is the set of wave-numbers for which S^^^ has -|-1 as 
an eigenvalue. This corresponds to a solution where no currents flow in the leads so 
that the conservation of current is satisfied on the internal bonds. 1 is in the spectrum 
of if 

Cg{k) = det [1 - S'^^\k)] = . (15) 
Eq. ()15|1 can be transformed in an alternative form 

. 1 r ,det[l -D(k)R] ^ ~ , r/'^Vf^ ^ ^ 

Ccfc =detl-p ^ ^^ = ; Rir-s i = Ri rs i + Sr s- TT- 16 

; [ ''Jdet[l-D(A:)i?] ' '"'l-pM ^ ^ 

which is satisfied once 

det[l - D{k)R] = 0. (17) 

In contrast to R, R is a. unitary matrix in the space of directed bonds, and therefore 
the spectrum is real. Eq. (|T7j) is the secular equation for the spectrum of the compact 
part of the graph [13 d ESj • 



4-2. Reaction-matrix formalism 

Eq. ()15|1 allows to compute the spectrum of the closed graph from the S"— matrix of 
the open graph. This relation can be inverted: it is possible to express the S"— matrix 
at an arbitrary wavenumber k in terms of all eigenvalues and eigenfunctions of the 
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20 



k 



Figure 1. A tetrahedron with a single attached scattering channel (left) is used to 
illustrate the reaction-matrix theory for graphs. The empty circle shows an auxiliary 
vertex with Neumann b. c. which is used for the derivation of Eg. 1201 On the right we 
compare the exact phase of the S-matrix of this graph (full line) to an approximation 
computed with the help of Eqs. fl^. (PT]l from the lowest 500 eigenstates of the graph 
(dots). 



corresponding closed graph. This is achieved by applying to quantum graphs methods 
which were developed in the context of nuclear physics: Wigner and Eisenbud jSU] 
introduced in the late forties the reaction matrix as a connecting element between a 
closed and the corresponding scattering system. The projection-operator formalism of 
Feshbach allows a systematic separation of the Hilbert space of a closed system 
from the embedding scattering system and on this basis Weidenmiiller and coworkers 
developed a theory expressing the S-matrix in terms of an effective non-Hermitian 
Hamiltonian [22] which contains the dynamics of the closed system and the coupling 
to the scattering channels. This effective Hamiltonian is the basis for the RMT results 
on resonance statistics that we will discuss in the following sections, and therefore it is 
important to understand how reaction-matrix theory and Weidenmiiller approach work 
in quantum graphs. 

In order to achieve this understanding we consider first the Green's function G of 
a closed graph. It can be defined by the infinite spectral sum 



where denote the nth eigenvalue and $„(a;b) the corresponding eigenfunction at 
distance Xh at the bond b. The graph will later be opened by attaching leads to the 
vertices 1 < j < M on which we assume Neumann boundary conditions (the b. c. at 
other vertices are arbitrary). We denote by Gjj/ the values of the Green's function on 
these vertices and introduce the Hermitian reaction matrix 



G{xb,Xb'] k) = Y, 



n=l 



oo 



(l>n{Xb)K{Xb') 
k^-kl 



(18) 




(19) 



We will show that K is related to the S-matrix via 



K 



i{i-sy\i + s) 



(20) 
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which is equivalent to 



S 



{I + iK)-\l -iK) . 



(21) 



For technical reasons it is easier to derive Eqs. ()2()j) . ()2H1 in a situation where the leads 
are attached to vertices j with valency Vj = 1 and Neumann b. c. An example is shown 
in Fig. Q As long as we consider Neumann b. c. this is no restriction of generality 
since the length 6 of the bond reaching a "dead end" at vertex j in the closed graph 
can be arbitrarily small, and it can easily be shown that in the limit 6 —>■ spectrum 
and eigenfunctions of the graphs with and without the attached arm coincide. Hence, 
after Eq. ()20|) is established with the help of auxiliary vertices, these vertices can be 
disregarded. Note that due to this construction p^^^ in Eqs. (fT2|) . (fT^ is zero: in the 
scattering system the vertices to which the leads are formally attached have perfect 
transmission and are effectively absent. 

There are two independent ways to arrive at Eqs. ()2()|1 . ()2H1 . The first simply repeats 
the corresponding arguments for billiard systems In fact the case of a quantum 

graph is completely equivalent up to the additional simplification that the attached leads 
do not allow for multiple transversal modes. However, here we will follow a different 
route which is more specific to quantum graphs as it relies on the exact representation 
of the corresponding S-matrix as a sum over classical trajectories (see Eq. ()14|l ). There 
is a similar expansion of the Green's function in terms of classical trajectories 



where t denotes a trajectory leading from x' to x while Bt, kit, and 9* are defined 
as in Eq. ()14j) . Skipping a more mathematical derivation of Eq. ()22j) we only point 
the analogy with the well known Green's function of a free particle Gq in one-dimension 
Go = (2iA;)~^e*'^'~^''. Eq. simply generalizes this expression by adding the additional 
paths from x' to x which are due to the multiply connected nature of the graph. 

Let us now consider the specific matrix elements Gjjr of the Green's function 
Eq. ()22|) . which appear in the definition of the reaction matrix Eq. ()19|). It is our goal 
to represent them in terms of the scattering trajectories contributing to the S-matrix 
element Sjf according to Eq. ((T3j). To this end we note that the contribution of each 
trajectory is the same in the two expressions up to an overall factor [2ik)~^ appearing 
in the Green's function. However, the set of contributing trajectories is different: in 
the open system we have p^^^ = meaning that no backscattered trajectories from the 
auxiliary vertices 1 < j < M contribute, while in the closed graph such trajectories 
do exist. In order to get their contributions also from Eq. (fT^ we consider multiple 
scattering at the graph with scattering trajectories reinjected an arbitrary number of 
times. Then it becomes indeed possible to represent all trajectories inside the closed 
graph as combinations of scattering trajectories and we find 



G{x, x'; k) 



1 



(22) 



2ik 



■X 




(23) 
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where n counts the multiple scattering at the graph. The first term (n = 0) corresponds 
to the zero- length trajectories from j to j. Note that Gjj' actually denotes the limit 
of G{x,x'; fc) as X — s> j and x' j' and that there are always two options to reach a 
point in the vicinity of a vertex j: with or without preceeding scattering at this vertex. 
Combining these options for j and j' gives the prefactor 4 (2 for n = 0) in Eq. (j2Sl)- 
After summing the geometric series in Eq. ()23|) we arrive at Eq. ((201) ■ 

Not only have we established in this way an analogue of the well known reaction- 
matrix formalism for quantum graphs, we have also given the first semiclassical 
interpretation of this theory: our derivation of Eq. ()20|) was entirely based on comparing 
the contributions of classical trajectories summing to an S-matrix element or Green's 
function, respectively. 

In Fig. H we demonstrate numerically, how Eq. (j2T|l can be used to compute the S- 
matrix of a scattering graph from its corresponding closed analogue. For the tetrahedron 
graph (see figure) we computed the lowest 500 eigenvalues and eigenf unctions and using 
Eqs. ()18j) . (PT|) we obtain an approximation for the 5'-matrix of the graph (see circles). 
For simplicity we consider the case M = 1 where the 5— matrix reduces to a single phase. 
The resulting data are in excellent agreement with the exact result (solid line) calculated 
from Eq. ()12|). We verified that the comparison is equally successful for multiple leads 
and other types of graphs. 

We continue by deriving the effective Hamiltonian for a quantum graph. To this 
end we rewrite Eq. (fTI?|) in the standard form 

oo T 

Kn'=^Y.W:,^^W^,,. (24) 

n=l '^n 

where 

Wn.Ak) = \fkhr'4'n,j (25) 

denotes the coupling of the internal level n to the scattering lead j. We see from Eq. (f^H|) 
that these coupling constants are proportional to the internal eigenfunction at the point 
where the lead is attached, as it should be expected. Following [2H1 we can now use 
some algebra to rewrite Eq. ()21|) in the form 

S = I -2-KiW^ —^—W , (26) 
E — rt 

where the S-matrix is expressed in terms of an effective Hamiltonian 7i, which has in 
the eigenbasis of the graph the matrix elements 

Unn' = kXn' ' tn{WW%^, . (27) 

Note that this Hamiltonian depends on energy via Eq. ()25|) . Whenever this dependence 
can be neglected, Eq. (j27j) can be diagonalized to find the resonance poles. In this regime, 
the predictions of RMT for the statistical properties of the resonances [211 12Z1 12H] are 
applicable. 
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Figure 2. Poles of the S*^^-* -matrix for regular Neumann graphs. The exact evaluated 
poles are indicated with (o) while (★) are the results of the perturbation theory H31|l : 
(a.) V = 5 and v = A and (b) V = 15 and v = 14. 



5. Statistical properties of resonance widths 

One of the basic concepts in the quantum theory of scattering are the resonances. They 
represent long-lived intermediate states to which bound states of a closed system are 
converted due to coupling to continua. On a formal level, resonances show up as poles 
of the scattering matrix S^^^^ occurring at complex wave-numbers Kn = kn — |r„, where 
kn and r„ are the position and the width of the resonances, respectively. From ()12|1 it 
follows that the resonances are the complex zeros of 

CgiK)=det{l-SBiK;A))=0 . (28) 

The eigenvalues of 5*^ = D{k; A)R are in the unit circle, and therefore the resonances 
appear in the lower half of the complex n plane. 

One important parameter that is associated with the statistical properties of 
resonance widths is the Ericson parameter defined through the scaled mean resonance 
width as: 

(7>. - (29) 

where {■)k denotes spectral averaging and A is the mean spacing between resonances. 
The Ericson parameter defines whether the resonances are overlapping ((7)^ > 1) or 
isolated {{■y)k < 1)- 

For relatively small values of the Ericson parameter a perturbation theory can be 
carried out in order to evaluate the resonance widths. Indeed, the difference 6R = R — R 
gets smaller as larger graphs are considered (for graphs with Neumann boundary 
conditions it is easy to see that the difference is of order That is, the leads are weakly 
coupled to the compact part of the graph. To lowest order, {6R = 0), the resonances 
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coincide with the spectrum of the compact graph. Let A;„ be in the spectrum. Hence, 
there exists a vector |?t,) which satisfies the equation 

D{kn)R \n) = l\n) . (30) 

To first order in 6R, the resonances acquire a width 



{n\ 


\D{K)5R\n) 


{n\ 


\L\ 


\n) 



(31) 



The above equation shows clearly that the formation of resonances is closely related to 
the internal dynamics inside the scattering region which is governed by Sb- 

To check the usefulness of Eq. (|3ip. we searched numerically for the true poles for 
a few scattering graphs and compared them with the approximation ()31|) . In Fig. |21we 
show the comparison for fully connected Neumann graphs with \^ = 5, 15 and A = 0. As 
expected the agreement between the exact poles and the perturbative results improves 
as V increases. 



5.1. Resonance width distribution 

An important feature of the distribution of the resonances in the complex plane can 
be deduced by studying the secular function Cg{.i^)- Consider Cg{i^ = 0)- If one of the 
eigenvalues of the matrix R (0) takes the value 1, Cg{k = 0) = and because of the 
quasi-periodicity of (g, its zeros reach any vicinity of the real axis infinitely many times. 
For f-regular Neumann graphs without magnetic vector potential {A = 0) the matrix R 
has indeed an eigenvalue 1, and therefore the distribution of resonance widths is finite 
in the vicinity of 7 = for these systems. For generic graphs, the spectrum of R is 
typically inside a circle of radius fimax < 1- This implies that the poles are excluded 
from a strip just under the real axis, whose width can be estimated by 

^ gap ^^i\f^max\^ / Lmax ; (32) 

where L^ax is the maximum bond length. The existence of a gap is an important feature 
of the resonance width distribution V{'y) for chaotic scattering systems. 

A similar argument was used recently in [21] in order to obtain an upper bound for 
the resonance widths. It is 

^max ~ ~2 In (I /imin I) /-^mm ; (33) 

where f^min and L^in are the minimum eigenvalue and bond length, respectively. 

The distribution of the complex poles for a generic fully connected graph with V = 5 
is shown in FigEti'- The vertical line which marks the region from which resonances are 
excluded was computed using (jHH). 

Random matrix theory can provide a general expression for the distribution of 
resonances. Specifically, Fyodorov and Sommers [2Z1 i2H] proved that the distribution of 
scaled resonance widths for the unitary random matrix ensemble, is given by 

r(M) ^ dj" \ (77r) i • ^ ' 
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Figure 3. (a) The 5000 resonances of a single realization of a complete V = 5 graph 
with A and M = V. The solid line marks the position of the gap jgap- (b) The 
distribution of resonance widths V{'j)- The solid line is the RMT prediction H34fl . The 
difference V^j) — VcuEil) is shown in the inset. 



where the parameter g = — 1 controls the degree of couphng with the channels 

(and it is assumed that Qi = g Wi = 1, ...M). 

In the limit of M ^ 1, Eq. (jH^ reduces to the following expression 



M t^^ M / _ / M 



for -7^<7< 



= ) 27r72 7r(3+l) ^ ' ^ ii{g-l) _ (^35^ 

I otherwise 

It shows that in the limit of large number of channels there exists a strip in the complex 
K— plane which is free of resonances. This is in agreement with previous findings 
inSl inni EZI- in the case of maximal coupling i.e. g = 1, the power law extends to 
infinity, leading to divergencies of the various moments of 7's. Using (jH^ we recover 
the well known Moldauer-Simonius relation |36l for the mean resonance width |2l 



\l)k = — • (36) 

The resonance width distribution for aV = 5 regular and generic graph is shown in 
Fig. Eb together with the RMT prediction, which reproduces the numerical distribution 
quite well. Fig. |3] shows a similar comparison for Neumann graphs. In this case, a 
clear deviation from the RMT prediction is observed. The main feature of V{'~f) is 
the relatively high abundance of resonances in the vicinity of the real axis which is 
of dynamical nature and conforms with the expectations (see related discussion above 
Eq. dSl). 

Although the absence of a gap in a strict sense is apparent in all the Neumann 
graphs that we have studied, still one can speak about a gap in a "probabilistic" sense. 
To show this, we have calculated the ratio of the resonance width standard deviation 
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Figure 4. (a) Resonance width distribution 7^(7) for fully connected Neumann graphs 
with {a.) A = and M = V = 5; {h) A ^ and M = V = 15. The solid hue is the 
RMT prediction H34() . In the inset we report the relative ratio -^y- vs. V for fully 
connected graphs with A = 0. As V increases, the distribution is concentrated more 
around its mean value and an effective gap is created. The dashed line has slope —0.25 
and is the best fit to the numerical data. 



A7 with respect to the mean resonance width (7)^ for various fully connected graphs. 
Our numerical data reported in the inset of Fig. |3] for various values of V suggest that 
^ ~ ^ with a ~ 0.25. 

For g ^ 1 (i.e. weak coupling regime) Eq. (jH^ reduces to the following expression 

where /3 = 1(2) for systems which respect (break) time-reversal symmetry, and r{x) is 
the Gamma-function. We notice that (j37p is the so-called distribution and can be 
derived independently using simple perturbation theory 

A simple graph that satisfy the weak coupling limit is the Neumann star graph. 
They consist of vq bonds and one lead, all of which emanate from a single common 
vertex labeled with the index i = 0. It is a simple matter to derive the scattering matrix 
g{M=i) ^]^[q}^ ig reduced to a single phase factor |2S1 

^(*' - - ;i::::lSit: p*^' 

In Fig. El we present our numerical results for V{'~f) for a star graph with vq = 20. 
The data are in excellent agreement with the RMT expectation given in Eq. ()38j) . We 
point out that in this case the coupling to the continuum is weak since (7 ~ 10 ^ 1 and 
therefore the x^-distribution with M = 1 is applicable. 



Statistical properties of resonance widths for open Quantum Graphs 14 

10^ 

P 10" 

10" 

10" 

10"' 

' ^ -? -1 n 1 

10 10 10 10 

Y 

Figure 5. The rescaled resonance width distribution 7^(7) for a star graph with 
vo = 20 (inset). The solid hne is the RMT prediction Eq. lpT7| . 




5.2. Resonance-trapping in graphs 

As mentioned at the end of Section|^2l the effective Hamiltonian Eq. is in general an 
infinite-dimensional and energy-dependent matrix and thus very hard to treat. However, 
there are situations where this Hamiltonian can be approximated with good accuracy 
by a finite and constant matrix. For example, this is the case, when the energy levels 
of the closed system are clustered, e. g. by an underlying shell or band structure. 

Let us calculate under this assumption the Ericson parameter Eq. (jSHI) of a group of 
N resonances of a graph at large /c 3> A. The resonances can be found by diagonalizing 
Ti. From Eq. ()27|1 we have for graphs with Neumann b. c. 

(7)^= -^I-trT^ 



n,] I 



= M^^^ . (39) 
A ^ ^ 

Note that diagonalization of Ti yields the resonances in the complex energy plane, while 

we are interested in the wavenumber k. This was accounted for in the first line of 

Eq. From the normalization of the eigenf unctions of the graph we deduce that 

the mean intensity of an eigenstate is the inverse of the total length C of the graph 

{\ipn,j\'^) = jC^^- Hence, if A was equal to the mean level spacing of all levels, tt/C, the 

Ericson parameter would be fixed to {'y)N ~ M/tt. However, we said above that the 

effective Hamiltonian can be used only if there is clustering of levels, meaning that the 

local mean spacing inside the cluster of N resonances is A ^ 7r/£. Therefore we have 

in such a situation not only (7) at ^ 1 but even 7 ^ M. 

Theoretical calculations, based on RMT modeling, predict that in this situation 

the statistical properties of the resonances change drastically: with increasing strength 
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Figure 6. Ring-shaped periodic graph with = 10 unit cells. 

of the coupling, M very broad resonances are formed, whereas the remaining N — M 
approach the real axis and become long hved (trapped) jHZl EH] • This counterintuitive 
phenomenon, which was recently demonstrated in experiments on microwave resonators 
[10] is called "resonance trapping" . 

Let us now illustrate this phenomenon for the case of quantum graphs. To this 
end we use a graph model consisting of identical unit cells forming a periodic ring as 
shown in Fig. El (we note that our numerical data were already included in the recent 
review by Dittes [32] )• Each unit cell consists of four vertices. Two of these vertices 
are T-shaped junctions with Neumann b. c. (A = 0). The other two are dead ends 
where we impose Dirichlet b. c. (A = oo). In one of the unit cells the dead ends have 
variable boundary conditions < Aieads < oo and are connected to a scattering channel. 
As long as the coupling to the channels is zero (A = cxd at all dead ends), the spectrum 
is organized in discrete bands, each containing N eigenvalues. These eigenvalues are 
seen in Fig. [3 as intersections of the solid lines with the real axis. We observe that 
the spacing within the bands can indeed be much smaller than the spacing between the 
bands. 

As soon as Aieads is finite, the eigenstates of the closed system are turned into 
resonances with a finite lifetime. In Fig. [7] we show the motion of resonances in the 
complex K— plane as we increase the coupling constant of the system by decreasing 
Aieads to zcro. One can clearly see the resulting redistribution of the S"— matrix poles. 
For particularly narrow bands, e. g. at ~ 4.2, we observe indeed two broad resonances 
with a width several orders of magnitude above that of the remaining poles, while for 
the broader bands at least one broad resonance can be distinguished. 

6. Conclusions 

Quantum graphs are used in our days a a paradigmatic model of quantum chaos. By 
connecting them with infinite leads one turns them into scattering systems. As a result, 
a scattering formalism can be developed, and an exact expression for the scattering 
matrix can be found. 

Due to the relative ease by which a large number of numerical data can be computed 
for the graph models, we have performed a detailed statistical analysis of resonance 
widths. A gap for the resonance widths has been obtained for "generic" graphs and 
its absence was explained for Neumann graphs. Finally, we have shown that the 
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Figure 7. Trajectories in the complex k— plane of the resonances of a ring graph with 
N = 20 unit cells as Aieads falls from oo to 0. For the maximum coupling to the two 
attached scattering channels (Aieads = 0), the location of the resonances is marked by 
squares. 



phenomenon of resonance trapping can appear also in quantum graphs. 
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